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Abstract 

We consider extensions of the standard model containing additional heavy 
particles ( "quirks" ) charged under a new unbroken non-abelian gauge group 
as well as the standard model. We assume that the quirk mass m is in 
the phenomenologically interesting range 100 GeV-TeV, and that the new 
gauge group gets strong at a scale A < m. In this case breaking of strings 
is exponentially suppressed, and quirk production results in strings that 
are long compared to A~^. The existence of these long stable strings leads 
to highly exotic events at colliders. For 100 eV < A < kcV the strings 
are macroscopic, giving rise to events with two separated quirk tracks with 
measurable curvature toward each other due to the string interaction. For 
keV < A < MeV the typical strings are mesoscopic: too small to resolve in 
the detector, but large compared to atomic scales. In this case, the bound 
state appears as a single particle, but its mass is the invariant mass of a 
quirk pair, which has an event-by-event distribution. For MeV ^ A ^ m, 
the strings are microscopic, and the quirks annihilate promptly within the 
detector. For colored quirks, this can lead to hadronic fireball events with 
~ 10^ hadrons with energy of order GeV emitted in conjunction with hard 
decay products from the final annihilation. 



1 Introduction 



The LHC has energized the particle physics community with the promise of new 
physics at the TeV scale. This is the scale where the origin of electroweak symmetry 
breaking and the solution of the hierarchy problem must lie. Most studies of physics 
beyond the standard model therefore involve minimal models directly motivated by 
these problems, most notably the MSSM. However, history teaches us that the true 
physics may be non-minimal, and the most striking experimental discoveries may not 
have any obvious connection to the "big questions." It is therefore important to look 
for any new physics that can manifest itself by the enhanced energy reach of a new 
accelerator, especially if it arises from a simple extension of the standard model. This 
is especially important at a hadron collider such as the LHC, where large backgrounds 
mean that finding a signal often requires knowing what to look for. 

The classic example of simple new physics not directly motivated by electroweak 
symmetry breaking is a Z' . This involves extending the standard model with a 
gauge group, plus a new Higgs sector that breaks the symmetry. The only 

parameter of the new Higgs sector relevant for phenomenology is the Z' mass, so the 
only parameters in the model are the coupling g' and mz' ■ In addition, there 

is a discrete choice of the charges of standard model fields under UiX)' ■ (We assume 
that some of these charges are nonzero, otherwise the quirks are not observable.) 
Although the Z' mass is not directly tied to electroweak symmetry breaking, the 
focus is on the phenomenologically interesting regime (very roughly mz' ~ TeV) that 
is not excluded by existing experiments, but may be probed at LHC. 

In this paper, we consider another equally minimal extension of the standard 
model. We assume that there is an additional unbroken SU(N) gauge group with 
some fermions Q, Q in the fundamental representation. (The qualitative features 
of the model are unchanged if the particles are scalars rather than fermions.) This 
model is parameterized by the mass of the new particles rriQ and the SU (N) gauge 
coupling, which can be parameterized by the scale A where it gets strong. In addition, 
there is a discrete choice of the standard model gauge quantum numbers of the new 
fermions. (We assume that some of these charges are nonzero, otherwise the Z' is not 
observable.) We assume that the mass of the fermions is in the phenomenologically 
interesting range (very roughly 100 GeV <, uiq < TeV) that is not excluded by 
existing experiments, but may be probed at LHC. 

New strong interactions have been considered often in particle physics, usually 
with strong interaction scales at or above a TeV. We instead consider the case where 
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A -C TeV, in particular 



A < ruQ. 



(1.1) 



We therefore refer to the new gauge interaction as "infracolor." Note that if Q is the 
hghtest particle in the fundamental representation of infracolor then it is automati- 
cally stable, since there is no lighter state with the same quantum numbers. We have 
learned recently that this model was first considered by L. B. Okun [1], who called 
the new particles "thctons." This model was considered as a limit of QCD in Ref. [2]. 
Ref. [3] also mentioned this model as an example of a "hidden valley" model. 

This paper will consider the phenomenology of these models with values of A 
ranging over many orders of magnitude (roughly 100 cV to 100 GeV). All these 
values are natural, since A is related to the fundamental gauge coupling go defined at 
a scale iiq by 



where b is the 1-loop coefficient of the SU (N) gauge coupling beta function. The scale 
A is exponentially sensitive to the value of g^, so each decade of energy is roughly 
equally likely. 

Cosmology places strong constraints on new light physics, even if it is weakly 
coupled to the standard model. However, if the reheat temperature is sufficiently 
low (T < GeV) the infracolor sector is never in thermal equilibrium, and there are 
no cosmological consequences. This shows that there are no model-independent con- 
straints from cosmology on this physics. If we assume thermal abundances for the 
new particles the cosmology is complicated, but may also be viable [4]. 

This paper will focus on the collider phenomenology of this model at the qual- 
itative level. This phenomenology of this simple model is surprisingly exotic. The 
reason is that breaking of the infracolor gauge string is exponentially suppressed due 
to the large Q mass. As we will see, this leads to very exotic phenomenology, so we 
call the new particles "quirks." ^ 

The coUider phenomenology of quirks depends crucially on the length of the 
strings. This is set by the scale where the quirk kinetic energy is converted to string 
potential energy. Since the typical event is not close to threshold, it has kinetic energy 
~ tuq and gives a string length scale 



This can also be motivated by the replacements "strong" —^ "string" , "quark" — > "quirk" . 



A = /xoe-'"'/'^o, 



(1.2) 




(1.3) 
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We will consider collider signals for string length scales ranging from the size of 
detectors (~ 10 m) to microscopic scales. 

This paper is organized as follows. In Section 2, we briefly discuss model-building 
issues such as naturalness and unification, as well as indirect constraints from precision 
electroweak data, cosmology, and astrophysics. In Section 3 we discuss production of 
quirks and strings. In Section 4 we discuss signals for macroscopic strings. In Section 
5 we consider annihilation of quirks catalyzed by the string. In Section 6 we discuss 
the signals of mesoscopic strings, those that are too small to be resolved in a detector 
but large compared to atomic scales. In Section 7 we discuss the collider signals from 
microscopic strings. Section 8 contains our conclusions. 



2 Models and Indirect Constraints 



In this section, we discuss model-building issues such as naturalness and unification, 
as well as indirect constraints from precision electroweak constraints and cosmology. 
This discussion is fairly standard, and our conclusion is that there are no strong 
model-independent constraints on quirks from these considerations. 



2.1 Coupling to the Infracolor Sector 

Because we assume that the scale of infracolor strong interactions is below the weak 
scale, the hadrons of the infracolor sector are kinematically accessible to existing ex- 
periments. However, the standard model is uncharged under infracolor, and therefore 
a quirk loop is required to couple the sectors. Since the quirks are heavy, this leads 
to highly suppressed couplings to the infracolor sector. 

The leading coupling between the standard model and the infracolor sector at 
low energies arises from the diagram of Fig. la. This gives rise to the dimension-8 
effective operator 



The 2-loop diagram of Fig lb can couple the infracolor gauge fields to dimension-3 
fermion bilinears, but these have an additional helicity suppression in addition to the 
additional loop suppression, and are therefore suppressed. For uiq > 100 GeV this 
operator is far weaker than the weak interactions, so production of infracolor gauge 
bosons at colliders with energy below the quirk mass is completely negligible. Probing 
this sector at colliders requires sufficient energy to produce quirks directly. 
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Fig. 1. Loop graphs contributing to the coupling of the standard model 
and infracolor sector. 



The operator Eq. (2.1) mediates infracolor glueball decay, for example to photons 
or gluons. The rate is of order 
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Note that this is very sensitive to both A and uiq. We have 



cr ~ 10 m 



7 VTeVy 



(2.2) 



(2.3) 



.50 GeV. 

We see that the infracolor glueballs can decay inside a particle detector for A > 
50 GeV, while the lifetime becomes longer than the age of the universe for A < 
50 MeV. 



2.2 Star Cooling 

Stars with temperature T > A can potentially cool due to emission of infracolor 
glueballs. Due to the rapid decoupling of infracolor interactions from standard model 
interactions in Eq. (2.1), we find that this does not give interesting bounds. 

We will focus on bounds from SN1987A, which has the highest temperature 
(T ~ 30 MeV) of the astrophysical systems used to constrain light particles. We 
can estimate the bounds by comparing to axion cooling, which constrains the ax- 
ion decay constant fa ^ 10^ GeV. For both the axion and infracolor, the dominant 
energy loss mechanism is nuclear bremmstrahlung. 

Below the QCD scale the coupling Eq. (2.1) gives rise to an effective coupling of 
infracolor gauge fields to nucleons: 



loTT^mn 



(2.4) 
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Here Aqcd ^ 1 GeV is the scale of strong QCD interactions. Factors of Att have been 
put in using "naive dimensional analysis" [5] . This is to be compared with the axion 
coupling 

jCeS ~ ^^aNN. (2.5) 

We therefore have 

rate of infracolor production / fl'^fl''^-^QCD V / 



QCD 



rate of axion production y IdTr^rriQ J / fa 



(2.6) 



which gives a bound of mg > 10 GeV. Although these estimates are very crude, the 
fact that the infracolor emission falls as l/mg means that the rate is highly suppressed 
in the interesting regime uiq ^ 100 GeV. 



2.3 Cosmology 

If the infracolor gauge interactions and/or the quirks have thermal abundances early 
in the universe, there arc stringent cosmological constraints. This paper will focus 
mainly on collider physics, so we make only some simple remarks here, leaving a more 
complete investigation to future work. 

The rapid decoupling of the infracolor interactions means that infracolor glueballs 
are not produced if the reheating temperature is sufficiently low. Assuming Trh ^ A, 
the condition for infracolor interactions to be out of equilibrium is 

^-(lISj)^^™^^. (2.7) 

which is satisfied for Trh ^ GeV. This is easily sufficient for nucleosynthesis at T ~ 
MeV, the highest temperature about which we have secure cosmological knowledge. 

This is not an entirely satisfactory solution to cosmology, since it requires dark 
matter and the baryon asymmetry to be produced at low termperatures. This is 
possible with e.g. McV dark matter [6] and low-scale baryogensis [7]. We can avoid 
exotic low-tcmpcraturc cosmology by having quirks decay to infracolored states that 
are sterile under the standard model. These decays can have lifetimes long compared 
to collider time scales, but short enough to avoid cosmological constraints. We will not 
discuss the details here. For the present discussion it is sufficient that low reheat tem- 
peratures are not in conflict with nucleosynthesis, so there is no model-independent 
constraint from cosmology. 
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2-4 Precision Electroweak Data 



Precision electroweak data constrains new physics at the TeV scale. However, if the 
quirks are in a vector-like representation of the standard model gauge group they 
can have a TeV mass term that does not break electroweak symmetry. Furthermore, 
virtual quirks are necessarily created in pairs, so there are no tree-level effects on elec- 
troweak observables. There is therefore no constraint on such models from precision 
electroweak data. 

2.5 Model Building 

Next we discuss the plausibility of this kind of new physics. The existence of addi- 
tional gauge groups with matter in bifundamental representations is a hallmark of 
brane constructions in string theory. As we will see the most natural quirk sector is 
vectorlike, which means that it requires no additional projections of the kind needed 
to obtain a chiral theory such as the standard model. A quirk/infracolor sector can 
therefore arise simply and naturally from string theory. 

In fact, in realistic supersymmetric theories there is already at least one set of 
vectorlike fields, namely the Higgs bosons. These must have a supersymmetric 
term" at the weak scale, otherwise we have either hght Higgsinos or no electroweak 
symmetry breaking. Any mechanism that generates the /x term can also generate 
a weak-scale mass for the quirks. This means that no additional assumptions are 
required to explain the origin of the quirk mass in supersymmetric theories. It is also 
trivial to preserve gauge coupling unification in supersymmetric theories by assuming 
that the quirks come in complete GUT representations. The simplest example is that 
the quirks are in a 

5 e 5 ^ (3, 1) 1 e (3, e (i, 2) i e (i, 2)_i . (2.8) 

In this model there is no tree-level Yukawa interaction that can split the masses of the 
doublet. These splittings will arise from loop graphs, and will be very small. There 
is also no tree-level interaction that allows either the color triplet or the electroweak 
doublet to decay to the other, so this model naturally has both colored and uncolored 
quirks. 

In fact, a quirk/infracolor sectors have already appeared in some model-building 
motivated by the hierarchy problem. Such a sector was proposed in Ref. [8] to give 
additional loop contributions to the physical Higgs mass in supersymmetry. Scalar 
quirks ("squirks") appear in models of "folded supersymmetry" [9]. 
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Small values of A are perfectly compatible with grand unification. As an example, 
we consider the MSSM with an SU{2) infracolor gauge group, with quirks in the 5® 5 
representation (see Eq. (2.8)). The infracolor beta function is equal to the color beta 
function at one loop, simple unification implies that the infracolor gauge coupling is 
equal to the QCD gauge coupling at the scale of superpartner masses. The scale of 
infracolor interactions is then of order 100 MeV. If the theory above the TeV scale 
has respectively 1, 2, 3 additional pairs of infracolor fundamentals, the infracolor scale 
is respectively MeV, 10 keV, 100 eV. 



3 Quirk Production and String Formation 



3.1 Absence of String Breaking 

The reason that infracolor gauge strings do not break was already discussed in the 
introduction. A virtual quirk-antiquirk pair costs energy at least 2mQ, and will have 
a typical separation of order thq^. This lowers the string potential energy only by an 
amount of order A'^/rriQ <^ 2mQ, so this process cannot go on shell (see Fig. 2). An on- 
shell transition requires eliminating a length of string of order AL ~ mq/ A? ^ A^^. 
The rate for this transition will be exponentially suppressed. 

This transition is closely analogous to the Schwinger mechanism of pair creation of 
charged particles by a weak electric field [10]. For charged particles with m » E^/"^, 
the rate for pair creation is 

V/V = ^e-™'/^. (3.1) 

Modeling a gauge string as a flux tube with area A and color electric field E, the 
string tension is cr ~ ii^^A, so we have 

r/L ~ -^e-^"^«/^. (3.2) 
For a string of length L ~ mq/a, we estimate E ~ ttA^ and obtain a lifetime 

r^^e-Q/^l (3.3) 
mq 

For h./mq < 10~^ this is already longer than the age of the universe for mq > 
100 GeV. 
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(a) 




Fig. 2. Schematic view of color flux for quirk separation for (a) r ^ 
and (b) r ^ A~^. String breaking (c) requires a quirk-antiquirk 
pair to be created, which costs energy 2mQ ^ A. 

3.2 Quirk Production 

Quirk production involves energy and momentum transfer of order mg ^ A and 
Aqcd, and is therefore a hard perturbative process. The total cross section for quirk 
production at Tevatron and LHC at leading order in perturbation theory are shown 
in Fig. 3. This does not include Sommerfeld enhancement due to attractive infracolor 
and/or QCD interactions [11]. This will increase the cross section near threshold, and 
need to be included in a more detailed study. The Coulomb interactions are familiar, 
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Fig. 3. Quirk production cross section at the Tevatron and LHC 



so we consider briefly the Sommerfeld enhancement due to the long-range infracolor 
interactions. These become relevant only when the string length is longer than A~^, 
which requires /? > (A/ttiq)^/^. The linear potential will be a large perturbation on 
the state if the potential energy changes signiflcantly in one de Broglie wavelength. 
We therefore compute the ratio of this change to the kinetic energy: 

^ . ;^!Z!!!2/3 . ^ i < M (34) 

We see that the effects of the long-range potential are always small enough to be 
treated as a perturbation, although they may be numerically signiflcant for the largest 
values of A. 

Returning to Fig. 3, we conclude that the cross sections are substantial up to 
several TeV for LHC. (Note that the cross section is proportional to A^ic-) Many quirk 
signatures are completely background- free (as we will see), so even a few reconstructed 
events may be sufficient for discovery. 
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3.3 String Formation 



The effect of non-perturbativc infracolor interactions and the formation of an infra- 
color string has many points in common with hadronization of heavy stable quarks 
in QCD, so we review this first. 

Imagine that there is a heavy (m ^ Aqcd) stable quark (or squark or gluino) in 
QCD in addition to the light quarks. Below the free quark threshold at 2m these 
can be produced in a Coulomb-like bound state (quarkonium). Formation of such a 
low- lying bound state requires that the quirk pair be produced just below the free 
threshold, i.e. 2m — £■ ~ Q;|(m)m. Wavefunction overlap factors give a suppression 
of the rate by additional powers of a^{m), and so the rate for the production of these 
bound states is much smaller than the production rate for unbound quarks. The quark 
production cross section is dominated by quarks with kinetic energy K = E — 2m ~ 
m >> A. In this regime, threshold effects are unimportant and the production process 
is perturbatively calculable in an expansion in 03(771). 

We now consider the hadronization of heavy stable quarks with ~ m. Kine- 
matically, it is possible that a large fraction of the kinetic energy is converted to 
light hadrons, resulting in a jet surrounding the heavy quarks. However, because the 
quark is very heavy its kinetic energy cannot be efficiently converted into production 
of light hadrons. The basic reason is that the strong interactions have a range of 
order Aq^j^, so once the heavy quarks are separated by a distance r ^ Aqqj^ the 
strong interactions become perturbative. It is traditional in heavy quark physics to 
refer to the non-perturbative QCD interactions as "brown muck" to emphasize how 
little we know about it. The size of the force exerted by the brown muck is of order 
AqcD) so the total energy transfered from quark kinetic energy into light hadrons is 
only of order 

A£; ~ FAr ~ Aqcd- (3.5) 

There is a tail at large that can be described in perturbative QCD by additional 
hard gluons. 

We now turn to quirks. The infracolor interactions effectively have infinite range 
bccaTisc of the infracolor string, and we might worry that the conversion of quirk 
kinetic energy to infracolor hadrons (glueballs) never stops. We consider events far 
from threshold {K ~ ""^q), for which the string length L ~ mq/K^ ^ A~^, long 
enough to be a well-defined object. In this case the string rapidly straightens out, 
approaching a configuration close to its local ground state. 

To understand this, it is helpful to restate in a somewhat formal way the obvious 
fact that well-separated QCD hadrons from heavy stable quark production do not 
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continue to lose energy to hadron emission. The point is that a state consisting of 
well-separated hadrons is locally (on scales of order Aqcd) a boost of the ground state. 
Let us apply this point of view to a rapidly-stretching infracolor string with heavy 
quirks at the ends. In the center of mass frame, the middle of the string has zero 
transverse velocity. A long QCD string is described by the Nambu-Goto action (see 
below), which has no longitudinal excitations. This string configuration is therefore 
identical to the ground state in the center of mass frame. Near the ends of the string, 
only the acceleration of the ends represents a departure from a boost of a ground 
state. The acceleration is given by 

a = ^ ^ « A. (3.6) 

Because the acceleration is very small on the scale A, there is no energy loss to 
infracolor radiation from the ends. The infracolor strings can be thought of as being 
close to the static limit mq — > oo. This is qualitatively different from the open strings 
of string theory, which have massless ends. 

The non-perturbative infracolor "brown muck" is therefore effective in radiating 
glueballs only when the quirk separation is of order A"^ or less. Similarly to the case 
of heavy stable quark production in QCD, this results in an energy of order A being 
radiated into infracolor glueballs during the production process. 



3.4 Dynamics of Quirks and Strings 

We now discuss the motion of the quirk-string system produced as described above. 
As long as we are considering excitations of the string with wavelengths long compared 
to A~^, we can use an effective description in which the string is elementary. This 
is analogous to the chiral Lagrangian describing pion interactions for energies small 
compared to Aqcd ~ GeV. 

Gauge strings are described at long distances by the Nambu-Goto action. This is 
not a priori obvious, since there are other universality classes of strings that break 
additional Lorentz symmetry. For a clear discussion of this point, see Ref. [12]. 
Strong numerical evidence that the long-wavelength fluctuations of the QCD string 
are described by the Nambu-Goto action was obtained in Ref. [13] . 

The action for a pair of heavy quarks connected by a gauge string can be written 

2 

S = -mg Y,JdTi-aJdA + S^^t, (3.7) 

1=1 

where dTi is the proper length of the worldline for quirk i, dA is the proper area 
element of the string worldsheet, and S'ext represents the effect of external forces. 
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Fig. 4. Definitions used in quirk equations of motion. 



Here a ~ is tlie string tension. In tlie variation witli respect to tlie quirk position, 
tliere is a surface term from the string action that generates the string force on the 
quirks. We therefore obtain the quirk equation of motion 



^9 - ^11 

1 — vf s -\ , v± 

- Vf 



+ ^ext, (3.8) 



where v is the quirk velocity, and v\\ and v± are the components of the quirk velocity 
parallel and perpendicular to the string: 

v\l = {v-s)s, v± = V — v\\, (3.9) 

where s is a unit vector along the string pointing outward at the endpoints (see 
Fig. 4). The second term in brackets is similar to a Lorentz force, and is required by 
relativistic invariance. 

The gauge string is a dynamical object with its own complicated equation of 
motion. However, if the quirks have no further interactions after they are produced 
{e.g. with matter in the detector) then in the center of mass frame the string remains 
straight. Therefore, the only long- wavelength excitations of the string arise from 
quirk interactions with matter. If the string force is much larger than matter forces 

i^ext<A', (3.10) 

then we expect that the string will remain approximately straight in the center of 
mass frame. The maximum force from either ionization or nuclear energy loss is of 
order (100 eV)^, so the straight-string approximation is guaranteed to hold only for 
A ^ 100 eV. Note that this translates to L -C 10 m, so all but the longest strings 
of interest in colliders can be approximated as straight. The full string dynamics is 
sufficiently complicated that it would be useful to check this by direct simulation. 

A potential concern is that interactions of the quirks with matter involve collisions 
with momentum transfer that may be larger than A. For relativistic quirks, the energy 
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and momentum transfer in these processes is of order 



Apion ^rue^ MeV, 

(3.11) 

Apnuc ~ Aqcd ~ GeV. 

We now ask whether this leads to the emission of infracolor glueballs. The important 
point is that only the quirk has electromagnetic or QCD interactions, so this energy 
and momentum transfer is to the quirk, not the infracolor string, which is sterile 
under the standard model. The change in the quirk velocity is of order 

- — < 1. (3.12) 

This is a small perturbation as seen by the infracolor interactions, and does not lead 
to the emission of an infracolor glueball. This is very clear if we consider a heavy 
stable quark, which is surrounded by QCD brown muck, but has no string attached. 
In this case, the perturbation is equivalent to the quark remaining at rest while 
the brown muck gets a velocity A^; in the opposite direction. This transfers energy 
AqcdAi'^ <S Aqcd to the brown muck. If this energy is smaller than the mass of the 
lightest hadron that can be emitted (a pion in this case), there is no transition and 
the process is elastic. For quirks, the total mass of a long string may be much larger 
than A, but glueball emission is a local process with a scale set by A~^. We therefore 
expect hadron emission to be suppressed as in the QCD case. 

4 Macroscopic Strings 

We now consider strings with lengths longer than the tracking resolution of a typical 
detector, very roughly L > mm. In this case, the quirk and the antiquirk appear 
as separated particles connected by a string. Strings much longer than a detector 
size will not have observable effects on the quirk trajectories, so we are considering 
mm ^ L ^ 10 m corresponding to 

100 eV ^ A < 10 keV (4.1) 

for niQ ~ TeV. 

4-1 Anomalous Tracks 

One obvious signature in this case is the anomalous quirk tracks in the case where 
one or both quirks are electrically charged. Because the string tends to accelerate the 
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Fig. 5. Anomalous tracks from quirks witli macroscopic strings. 

quirks toward each other, we can have events such as those depicted schematically 
in Fig. 5. In these events, the curvature of the tracks is qualitatively different from 
the curved track of a particle in the magnetic field of the detector. For example, a 
magnetic field along the beam direction curves tracks only in the r-0 plane, while quirk 
tracks generally have curvature in the r-z plane. Therefore, unambiguous observation 
of only a single event of this type is sufficient for discovery of macroscopic strings! 

Do quirks annihilate when the string force brings them back together? For the 
case of macroscopic strings considered here, this is highly suppressed by the fact 
that annihilation requires the quirk to be in a state of relative angular momentum 
£ ~ 1, while interactions with matter change the angular momentum by much larger 
amounts due to the long lever arm. Even a single ionization interaction gives 

,4.) 

The infracolor "brown muck" surrounding the quirk has a much larger cross section 
of order A~^, and can therefore interact for angular momenta ^ < mg/A. A single 
ionization interaction changes the angular momentum more than this for A < MeV. 
We conclude that quirks with macroscopic strings do not annihilate. 

The difficulty in detecting quirks with macroscopic strings is that triggers and 
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track reconstruction algorithms are designed for conventional tracks, and will likely 
miss these events altogether. Defining an efficient trigger for these events that has 
low background from standard physics and instrumental noise is worth further inves- 
tigation. A simpler strategy is to focus on events where the quirk pair is produced 
in association with one or more hard jets or photons (see Fig. 5b and 5c). Standard 
reconstruction algorithms will fail to reconstruct the quirk tracks, resulting in missing 
Pt balanced by jets or photons. If such events are discovered, careful examination 
of the signal events in the missing pt direction can reveal the presence of "quirky" 
tracks. 



Do quirks stop in the detector? The stopping of heavy stable charged and/or strongly- 
interacting particles has been extensively studied [14], with the conclusion that typi- 
cally a significant fraction do indeed stop inside the detector. For quirks there is an 
additional complication from the string interaction. In order for quirks to come to 
rest, they must become bound to the lattice in the detector material. If the string 
force is stronger than the forces that bind the quirks to the lattice, they will continue 
to be dragged by the string. 

We first consider possible final states of quirks bound to the lattice. The binding 
mechanism depends on the standard model quantum numbers of the quirks. If quirks 
are electrically charged but uncolored, they can be electronically bound to the lattice 
similarly to ordinary nuclei. This is particularly clear for positively charged quirks, 
which can share a lattice site with an ordinary nucleus since there is no constraint 
from the exclusion principle. Negatively charged quirks will experience an electrical 
potential with opposite sign, and it is reasonable to assume that they will also find 
a stable local minimum. If quirks are colored, they will form quirk hadrons whose 
charge may change with time because of inelastic strong interactions that change the 
valence quark structure. We expect quirk hadrons to bind efficiently with nuclei, and 
these can also become stuck in the lattice. 

In all of these cases, the binding energy of the quirk (or quirk-nucleus bound state) 
to the lattice is of order eV, and the typical size of a potential well is of order A. 
(The binding energy is set by the electron mass, and is independent of the mass of 
the heavy particle. For example, the binding energies for heavy and light nuclei are 
all several eV.) Therefore, the force required to remove a bound quirk from the lattice 
is of order 



4-2 Stopping Quirks 




eV 



(100 eV)l 
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If the string force is larger than this, the lattice cannot bind the quirk and it will not 
stop. 

Even if A <^ 100 eV, the string force gives the quirks substantial kinetic energy, 
making it more unlikely for them to bind to the lattice. The binding energy of a 
nucleus in the lattice is of order eV, and it is reasonable to assume that the lattice 
cannot absorb energy larger than this without breaking. Therefore, a quirk nucleus 
will not bind with the lattice if its kinetic energy is large compared to eV. This 
requires 

/^.lO-(^)". (4.4, 

For such small values of f3 ionization forces are described by the theory of Fermi and 
Teller [15], extended by Lindhard [16]. We have 

Fion ~ A-lP, Ao ~ keV. (4.5) 

The ionization force for different nuclei in the same material vary over about an order 
of magnitude, suggesting an uncertainty of an order of magnitude in Aq. Balancing 
this against the string force gives a terminal speed 

Imposing Eq. (4.4) then gives 

A.ev(^)-. (4.) 

This bound is proportional to Aq, so there is an uncertainty of an order of magnitude 
in this estimate. Despite this uncertainty, it seems unlikely that quirks stop in the 
detector even for the smallest values of A of interest. 

If one quirk stops in the detector, the other will eventually lose its kinetic energy 
and annihilate with it. If both quirks stop, there is a string stretched between them. 
This string can interact with strings of subsequently produced quirks, producing 
even more bizarre events. One can also imagine releasing such quirks by e.g. melting 
the material in which they are trapped, and looking for the subsequent annihilation. 
These are amusing possibilities that might be worth taking seriously if more a detailed 
study indicates that large numbers of quirks in fact stop in the detector. 

4.3 Quirk Annihilation 

Although the string force tends to prevent the quirks from stopping, it also tends 
to bring them together to allow them to annihilate. In particular, sufficiently slow 
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charged quirks (or quirk nuclei) will reach the terminal speed given by Eq. (4.6). 
The subsequent motion of the quirks is damped, so these quirks can find each other 
more efficiently. Of course, this mechanism only works if both quirks are electrically 
charged. 

The linear form of the damping force holds up to velocities of order a. For 
such small velocities slow-moving colored quirks can bind with nuclei, making them 

effectively charged and subject to the mechanism considered here. The maximum 
ionization force is of order (100 eV)^, so this mechanism only works for A < 100 eV. 
Ionization energy loss is a good description as long as the separation of the quirks is 
larger than atomic distances of order A. Even for distances smaller than Athe energy 
loss is more complicated, but we expect charged quirks to exchange energy efficiently 
with electrons. We therefore assume that they annihilate on a time scale relevant for 
collider searches. 

An interesting question is the distribution of annihilation events in the detector. 
We have made a crude simulation of this using the straight string approximation. We 
include ionization energy loss as a continuous force, approximating the detector as 
solid iron. We also include a crude approximation to nuclear energy loss, although 
that does not really affect our results. Wc assume that all quirks that reach terminal 
speed and come close together in the detector annihilate sufficiently rapidly to be seen. 
An example of our results arc shown in Fig. 6. Note that most of the annihilations 
take place near the beam. This is easy to understand. Most of the events where both 
quirks become damped arise from events where quirks are produced nearly back-to- 
back in the central region of the detector. In such events the quirks will have speed 
less than a at the turning points, and therefore they become damped there. Their 
subsequent motion is essentially constant velocity toward each other, and they meet 
near the beam axis. 

The distribution of annihilation events is very different that of late-decaying par- 
ticles stopped in the detector [14]. More realistic simulations should be done to check 
the distribution of these events. Another difference is that most examples of late- 
decaying particles that have been discussed in the literature decay partly to missing 
energy while quirks will annihilate to visible energy in most modes. 

Another aspect of quirk annihilation that can in principle give a signal is the 
ionization track of the damped quirks before they re-annihilate. The ionization is 
large compared to typical particles, but the track is very slow (/3 ^ a). Presumably, 
it will therefore generate "stub" tracks in many events that are triggered for other 
reasons, and these stubs can in principle be connected. Since these tracks lead to 
annihilation events (assuming that the timescale for energy loss is sufficiently short) 
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m =800 GeV A = 90 eV at LHC 




z(m) 



Fig. 6. Results from a crude simulation of position of quirk annihilation 
events relative to collision point. 



so one can start looking for them there. 

The previous discussion assumes that both quirks are electrically charged, so that 
they both experience ionization forces. If the quirks are colored, their charge state 
may change on a distance scale given by the nuclear mean free path (~ 10 cm in iron) 
complicating the phenomenology further. One other case that bears mentioning is the 
case of uncolored quirks where one is charged and the other is neutral, e.g. produced 
by s-channel W exchange. In this case, the charged quirk can become damped, while 
the neutral quirk will not interact with matter. In this case, the charged quirk will 
be driven by the invisible neutral quirk which can have a much larger amplitude of 
motion. This can result in truly bizarre charged tracks such as the one illustrated in 
Fig. 7. Since the damped quirk is moving very slowly (/3 <^ a) these events will be 
very difficult to detect. 
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Fig. 7. Highly exotic track resulting from an electrically charged quirk 
(solid track) becoming damped in the detector, while its neutral partner 
(dashed track) drives its motion. The neutral track will be unobserv- 
able. 



5 Prompt annihilation 

We now consider in more detail the question of quirk annihilation, which is very 
important for the phenomenology of microscopic strings. The momentum transfer 
in the annihilation process is of order mg, which means that the quirks must come 
within a distance of order in order to annihilate. Equivalently, the cross section is 
dominated by partial waves with relative angular momentum I < p/mg < 1. Because 
the maximum quirk separation L is much larger than the microscopic scales rriQ^ and 
Aqqj-,, there is a large lever arm with which interactions with matter can change the 
angular momentum. However, if the string is sufficiently short matter effects are not 
important (we will be more precise about this below). In this section we analyze 
annihilation of quirks in the absence of matter effects. 

A crucial question is the rate of transfer of energy and angular momentum from 
the bound state. An important feature is interactions of the non-perturbative "brown 
muck" surrounding the quirks, from infracolor and/or QCD interactions. The cross 
section for these interactions is much larger than the hard annihilation of quirks, and 
may change the energy and angular momentum of the system, thereby suppressing 
annihilation. We also consider the effects of radiation as a mechanism of losing 
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energy and angular momentum, and we argue that this is generally unimportant as 
a mechanism of energy loss. 

5. 1 Quirkonium 

Quirk pairs produced near threshold can form a Coulomb-like "quirkonium" bound 
state. Formation of such a low-lying bound state requires that the quirk pair be 
produced near threshold, i.e. \E — 2mQ\ < a'^niQ, where a is the infracolor gauge 
coupling, or the QCD gauge coupling if the quirks are colored and A < Aqcd- These 
bound states will annihilate promptly into pairs of gluons or quarks (for colorful 
quirks). This signal has been considered for stable gluinos in Ref. [17] where it 
was found to be less sensitive than searches for unbound gluinos. We expect the 
result for quirks to be qualitatively similar, in that signals for highly excited quirks 
{E — 2m.Q ^ a^mg) will be more sensitive. 

5.2 Highly Excited Bound States 

We are interested in the majority of events that produce quirks that are not close 
to threshold, i.e. E — 2mQ ~ toq. As discussed in Section 3 quirk production is 
essentially perturbative, so quirk pairs are produced in a state of relative angular 
momentum ^ ~ 1. The subsequent infracolor (and possibly QCD) "hadronization" 
stage does not strongly affect the energy and angular momentum of the quirks, so the 
quirk pair still has £ ~ 1 even when the quirks have large separation {e.g. r ^ A~^). 

If there are no further interactions that change the quirk angular momentum, then 
the quirk pair will not have a well-defined angular direction even if it is macroscopic! 
If the quirks interact with matter {e.g. in the detector) their angular position will 
certainly be determined, but if we are interested in cases with sufficiently large A then 
the typical length of a quirk string will be small {e.g. < A) and matter interactions 
cannot "measure" the angular position. In this case, a collider will create the quirk 
pair in a "Schrodinger cat" state with large-scale (r ~ mg/A^ ^ A"-*^) quantum 
correlations. 

Such a situation is not familiar in particle physics, and we will proceed cautiously. 

In the end, many of the results can be understood from a simple classical picture, 
but we will derive the results using WKB wavefunctions to take into account the 
important quantum-mechanical aspects of these states. 
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5.3 Wavefunction Overlap 



Wc now consider the probability that a highly excited quirk pair can be found suf- 
ficiently close together to re-annihilate. This is a standard problem in quantum me- 
chanics, and we review it here to set the stage for the subsequent discussion. Highly 
excited states can be described using the WKB approximation. The quirk annihi- 
lation probability is proportional to the probability to find the quirk pairs within a 
distance of order rng^ of each other. We will estimate this probabihty using non- 
relativistic quantum mechanics and simple approximations that are sufficient for our 
purposes. 

We begin with the case £ — 0. We denote the radial Schrodinger wavefunction by 
■0(r) and define the reduced wavefunction by 

y(r) = (6.1) 



/47rr 

This satisfies the normalization condition 

i^/jyWI^ (5.2) 

Jo 

and the boundary condition 

1/(0) = 0. (5.3) 
The time-independent Schrodinger equation can then be written 

y"(r) = -k\r), (5.4) 

where 

k{r) = ^^^X - V{r) . (5.5) 

Here K = E — 2/xq is the kinetic energy, and /xq = mQ/2 is the invariant mass of 
the reduced system. We temporarily keep % ^ 1 to keep track of the classical hmit. 
We are interested in the case of a linear potential, but we will see that the important 
results of this section are independent of the details of the potential, so we will keep 
it general. 

We approximate the wavefunction of this system by 



y(r) ~ , sm 



f dr'kir') 
Jo 



0{r^^-r). (5.6) 
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where rmax is the classical turning point, i.e. 

Hr^s.) = 0. (5.7) 

For r ^ Tmax this is the WKB wavefunction, and has the correct boundary condition 
at r = 0. The boundary condition at the classical turning point is only crudely 
approximated, but this will not strongly affect the probability to find the particle 
near the origin. In this approximation, we can compute the normalization constant 
as 



1 = ICP / dr — sin^ / dr'k(r' 
Jo k{r) iJQ 



(5.8) 



For highly-excited states, we are averaging over many periods with a slowly- varying 
potential, so we can replace sin^ by its average value \. This gives 

1 = ^ , . (5.9) 

-^0 ^JK - V{r) 

This is directly related to the time for a classical trajectory to go from r = to 



/■''max f]r 

/ , (5.10) 

Jo JK-V{r) 



For a hnear potential 



V{r) = ar (5.11) 



we have 

T=^. (5.12) 

(T 

However, it is more insightful to leave the results in terms of T as we will see. We 
therefore have 

i^r = |^- (5.13) 

We now estimate the probability to find the quirks within a distance Tq of each 
other: 

Prob(r < ro) = / dr \y(r)\^. (5.14) 
Jo 
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Near the origin, the wavefunction is oscillating with the de Broglie wavelength 

A„ = (5.15) 

For To ^ Ao the integral averages over many periods, and we can again replace the 
sin^ term by its average value of \: 

Prob(r<ro).i|Crf J^ = ^. (5.16) 



where 

o / ; (5.17) 

2 Jo - V{r) 



At 




is the classical time to go from r = to r = Tq and we have used Eq. (5.13) to 
eliminate C. This is the result familiar from quantum mechanics textbooks that in a 
highly excited state the probability to find a particle at the origin is proportional to 
the fraction of time that a classical orbit spends there. 

In the opposite limit Tq <S Aq we can use the approximation 

sin / dr' kir') ~ sin /cq'^ — ^o^, (5.18) 
Jo J 

where Uq = 2t:/\q. Since the wavefunction at the origin is 

m = (5.19) 

this gives 

Prob(r < ro) = \C? / kor^ = 1-0(0)^1^, (5.20) 

where V — ^'kTq is the volume of the region of interest. This is the result familiar 
from positronium and quarkonium physics that the £ = annihilation probability is 
proportional to the wavefunction at the origin. 

The results Eqs. (5.20) and (5.16) are very different parametrically. For Tq ^ Aq 
the result is classical and therefore independent of %, which is not the case for tq <^ Aq. 
Also, for ro S> Aq the probability goes as ro (since At ~ ro/?; where v is the velocity 
of the classical trajectory near the origin), while for ro ^ Ao the probability goes as 

'o- 
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We will be mainly interested in the limit K <^ jdQ, where Tq -C Aq. From 
Eqs. (5.20) and (5.13) we have 



|,,,0)|^ = ^ = ^, (6.21) 

I'^v J\ 47rT' ^ ^ 

where we have expressed the result in terms of the physical quirk mass and the velocity 
of a single quirk at production, given by 

P-M". (5.22) 

For the majority of events /3 ~ 1 and therefore |^(0)p ~ T~^. As we will see, 
this means that for highly excited quirk bound states there is a definite annihilation 
probability per classical crossing. Eq. (5.12) shows that T oc /3, so the probability is 
nonzero at threshold. 

It is straightforward to include the effects of nonzero orbital angular momentum. 
The angular momentum barrier means that classically the particles have a distance 
of closest approach given by 

e 



^2mQK 



£Xo, (5.23) 



where i is the angular momentum. (We are setting h = 1 again.) For r < Vq <^ Xq < 
fmin we can use the approximation 

yt{r) ~ Cer'+\ (5.24) 

We can determine the coefficients Q by matching onto the wavefunction for r ^ r^in- 
Since rmm ^ Xq, the sine function in the wavefunction is of order 1, and we have 

Viirn.^) - - Qr'^^\ (5.25) 

With this approximation we obtain 

P„b(. < ,) . _L_ (^) (,26) 

For To ~ m'o' this is suppressed compared to the i — case by 



Prob,^o(r < mg^) 1 ( pV^' 
Prob,=o(r < mQ^) ^ £ \£ ) 



(5.27) 



This suppression means that annihilation is dominated by small £. 
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5-4 Annihilation Rates 



We now use the results above to compute the quirk annihilation rates. For now 
we neglect the effects of interactions with matter, non-perturbative interactions, and 
radiation. We work in the highly excited regime a^niQ <^ K <^ mq. Our results 
should be approximately valid for K < ttlq, the regime where the majority of quirk 
pairs are produced. In this regime, the state is sufficiently excited to use the WKB 
approximation of the previous subsection, but the de Broglie wavelength of the quirk 
is larger than the distance Tq ~ m^^ over which the annihilation takes place. As we 
reviewed above, the probability to find the quirks near the origin is dominated by the 
£ — partial wave, and is proportional to |'0(O)p. The density of particles within the 
range of the annihilation cross section is therefore 1-0(0) p, and the annihilation rate 
is 



^(JVrel, (5.28) 

where a is the annihilation cross section and v^ei is the relative velocity of the quirks 
in their center of mass frame. The wavef unction at the origin is given by Eq. (5.21), 
and is proportional to l/T, where T is the classical time for the quirks to go from 
r = to r = Tmax- The annihilation probability per classical crossing is therefore 

P = 2TT= "^avrei. (5.29) 

ZTT 

Note that all dependence on the potential has dropped out, so the time scale for 
annihilation is set by the classical crossing time. 

The annihilation cross sections can be computed perturbatively. To get numerical 
factors right, note that the spin of the quirks is not correlated, and so we must average 
over intial spins. Similarly, the QCD color of colored quirk pairs is uncorrelated, so 
we average over quirk colors. On the other hand, quirks are in an infracolor singlet 
state because they are connected by an infracolor string. 

For quirks carrying QCD color, we then have 

r«Q-9») = ^5^W0)P. (5.30) 

Z ( I tin 



V{QQ^uu)^^-^'^\m\\ (5.31) 

where crVrei has been replaced by its threshold value. The QED annihilation processes 
are 

r(QC-77) = ^5^W(0)P, (6.32) 
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WQ ^ 7* ^ e+e-) = —1^(0)1^ (5.33) 



where eg is the electric charge of the quirk and Nq is the number of QCD colors of the 
quirk (so Nc = 1 if the quirks are color singlets). We have neglected the contribution 
from Z boson exchange, which gives a small correction. There are similar expressions 
for annihilation through a 1^ in the case where the electric charge of the quirks differs 
by one unit. 

For quirks carrying both QCD color and electric charge there is a potentially 
interesting mixed annihilation to gluons and photons with rate 



r«Q-97) = ^^'^W0)p. (5.34) 



This motivates searches for photon-jet resonances at colliders. 
There is also annihilation to infracolor gluons, which gives 

r{QQ ^ infracolor) = ^1^^(0)1^. (5.35) 

For A < 10 GeV, the infracolor glueballs are stable on collider scales and this is an 
invisible decay. The gauge couplings in Eqs. (5.30)-(5.35) are to be evaluated at a 
renormalization scale mg. For the infracolor coupling, we approximate 

Here we estimate the scale where the perturbative coupling blows up as Aic/4. (In 
QCD, this scale is ~ 250 MeV, while the scale of strong interactions is ~ GeV.) 

The annihilation probability per classical crossing time is important in comparing 
the annihilation rate with other energy loss mechanisms. For colored quirks, the 
annihilation into quarks and gluons dominates. The probability of annihilation per 
classical crossing is 

Pqcd = 2TrQCD = ^ «3/3 ~ ^, (5.37) 

where we have assumed that the quirks can annihilate into all 6 quark flavors, and 
we assumed A^ic = 3 and i3 ~ 1. For uncolored quirks with opposite charge, the 
probability per classical crossing to annihilate is 

Pqed = iVic(6e| + e^a^^ ~ ^, (5.38) 
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where we have used cq = 1, Nic = 3 for the numerical estimate. The annihilation 
probability per classical crossing to annihilate to infracolor gluons is 

^'■^^ ^ ^A/i^'^'c^ ^ 160 [l-0.121n(Aic/GeV)]2' ^^"^^^ 

where we assume uncolored quirks with A^ic = 3 and mq ~ TeV. We see that the 
annihilation into visible final states is significant even for uncolored states with large 
values of A. These results will be useful in assessing the probability that quirks 
undergo prompt annihilation. 

Another important quantity is the branching ratio for colored quirks to annihilate 
to leptons and photons. We have 



2 ^,2 



TjQQ ^ 9 ej^a 

TiQQ^jj) 68 a 



4 „,2 



1.4 X 10-^ (5.40) 



TjQQ ^ 11) _ 9 4^ .,3,, 10-5 .54l^ 

where we have taken eg — ^ for the numerical values. These branching ratios are 
discouragingly small.^ On the other hand, the branching ratio to photon plus jets is 



r(QQ ^gi) _ 6 



3 X 10"^ (5.42) 



r(gg^jj) 17 as 

where we again take cq = |. This is somewhat more encouraging. 

The accuracy of the estimates above can easily be improved by incorporating the 
behavior of the cross section at threshold and relativistic effects. We leave this to 
future work. 

The annihilation probabilities computed above are relevant when there are no 
interactions that can change the angular momentum of the quirks. We now consider 
these interactions to see whether they are in fact negligible. 



5.5 Non-perturbative QCD Interactions 

We now consider the effects of non-perturbative QCD interactions on colored quirk 
annihilation. Colored quirks are surrounded by a cloud of non-perturbative QCD 
"brown muck" with size i^had ~ ^qcd- Interactions between the brown muck of 

^These branching fractions arc lower than the corresponding ones for the Upsilon decays mainly 
because the initial state is not a color singlet. This opens additional colored channels and enhances 
the strong decay rate. 
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the quirks are important because they have a larger cross section than the hard 
annihilation processes considered previously. Although they do not result in the 
annihilation of the quirk pair, they can change the angular momentum of the quirk 
pair and affect the probability for hard annihilation. 

We assume that L ^ -Rhad, so that we can treat the quirk hadrons as well- 
separated particles moving under the influence of an infracolor string. This requires 

E Z T . (5.43) 

-'4-QCD 

For example, for A ~ GcV this requires only E > GeV, a mild requirement for hard 
production. For A > 10 GeV, this starts to be a significant constraint, and our results 
will be qualitatively reliable at best. 

The brown muck will interact only when the quirks come within a distance of 
order ii!had- We are interested in inelastic processes {e.g. pion emission) that can 
change the energy and angular momentum of the bound state. The typical energy 
transfer from the bound state can be estimated from 

AE - FAr - A^cD^had ~ Aqcd- (5.44) 
The momentum transfer can in principle be larger if the quirks are moving slowly: 

Ap ~ FAt - A|cD— ~ — ■ (5.45) 



However, for pion emission Ap ~ AE, so energy and momentum transfer are of order 
Aqcd- We expect such processes to have a geometric cross section of order 7rit!^ad) 
since there is no small parameter suppressing the interaction probability.^ Note that 
the large kinetic energy carried by the quirk is not transfered in the interaction, and 
does not suppress the interaction probability. 

A geometrical cross section is equivalent to saturating unitarity for all partial 
waves up to 

^max ~ rnqvRh^d ~ (mQ£;)^/^i?had- (5-46) 

Unless we are very close to threshold we have fmax ^ 1, and so the partial wave cross 
section will approximately saturate unitarity unless the angular momentum is very 
large. This means that the interaction will take place with of order unit probability 
whenever the quirks come within a distance of order i?had or less. We therefore expect 



"^If the number of QCD colors is regarded as a large parameter, the probability for interaction is 
of order I/Nq- We will neglect large A^c effects in the following. 
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that a brown muck interaction transfering energy and momentum of order Aqcd will 
take place roughly once every crossing time, as long as £ < £max- 

We can understand this result using a simple quantum-mechanical model. We 
model the brown muck as a particle (a constituent quark) of mass ~ Aqcd bound to 
each quirk by a potential that represents the effects of the QCD interactions. The 
wavefunction for the system is then a function of the relative coordinate of the quirks 
f and the coordinates of the constituent quarks relative to the associated quirk pi 2- 
The wavefunction is assumed to take the approximate form 

^(r, Pi, P2) ~ ^(r)x(pi)x(p2), (5.47) 

where '4}{r) is the quirk wavefunction, and x(p) is the consituent quark wavefunction. 
This factorized form is justified for r ^ -Rhad where the quirks are well separated. The 
constituent quark wavefunction x(p) is nonzero only for p < i?had- It is convenient to 
normalize it so that 



d'p\x{p)\'^Ri^. (5.48) 



We then have (see Subsection 5.3) 



w(r) ~ . sm 



f dr'kir') 
Jo 



(5.49) 



We are using an £ = wavefunction, which will have the right qualitative behavior 
as long as i ^ ^max- Normalizing the wavefunction gives 

\Cf ~ jr^, (5.50) 

''■-"■had-' 

where T is the classical crossing time. Here we make the same approximations as 
previously for the quirk wavefunction. We can then compute the probability that the 
brown muck particles are within a distance Rhad from each other: 

Prob(|f+pi-p2|<i?had) = ^ jd'r^ jd^p^\x{pi)\^ jd%\x{p2)? 

x^(i?had-|r + pi-p2|) (5.51) 

where v ~ {E/rriQ)^^'^ is the classical quirk velocity at the origin and we have again 
used the fact that the Compton wavelength of the the heavy quirk is much smaller 
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than -Rhad- This is the fraction of the time that the quirks are within a distance i?had- 
To find the reaction rate, we must find the density of incident particles over the range 
of the interaction. This is 

SO the reaction rate is 

pvar^ vRl^ ~ ^. (5.54) 

We are again led to the conclusion that these interactions occur roughly once per 
classical crossing. 

One important effect of these interactions is that it changes the angular momentum 
state of the quirk pair. The quirks are produced in a state with angular momentm 
^ ~ 1, ie. a highly spherical quantum state in which the angular position of the quirks 
has nearly maximal uncertainty. The hadrons that are emitted eventually interact 
with matter far from the detector, and therefore can be thought of as having a definite 
direction. The fact that the angular momentum state of the hadrons is entangled with 
that of the quirk pair means that this reduces the quantum uncertainty in the the 
angular direction of the quirk pair. In the traditional textbook language of quantum 
mechanics, the angular position of the quirks gradually becomes "measured" by the 
repeated "measurement" of the pion angular positions.^ A proper treatment of this 
process using the ideas of quantum decoherence is beyond the scope of the present 
work, and we will only make some simple estimates here. 

The angular momentum transfered to the emitted hadrons in a single brown muck 
interaction is of order 

Ai ~ i?hadAp ~ 1. (5.55) 

Assuming that the interaction is equally likely to raise or lower the angular momen- 
tum, we have i ~ after N such interactions.^ The hard annihilation cross section 
falls rapidly for i > few, so there is a competition between hard annihilation, which 
wants to eliminate the bound state in a small number of classical crossing times, and 
the non-perturbative QCD interactions, which tend to increase the average angular 
momentum, and therefore suppress hard annihilation. 

"^We are neglecting possible interactions of the quirks with matter, which would directly "mea- 
sure" the angular position of the quirks. This discussion is therefore applicable to the case where 
the string is sufRciently short that matter interactions are unimportant. 

^The angular momentum is positive semi-definite, but this is taken care of in the random walk 
by simply identifying ±£. 
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We can illustrate these points with a simple quantum mechanical toy model. We 
work in 2 dimensions, where the angular momentum eigenstates are simply e""^, 
where 9 is the polar angle and m is an integer. We can simplify the model further 
by restricting the particles to a circle, so there is no radial wavef unction to worry 
about. We assume that there is a process by which a "quirk bound state" in an 
angular momentum m state emits a "pion" that also lives on the circle. The 1-particle 
wavefunction therefore makes a transition to a 2-particle wavefunction 



ooe 



+ a.e'^^^'^'e-''' + a^.e'^^-^^'e''' + • • • . (5.56) 



Here 9' is the angular coordinate of the emitted pion. The transition conserves angular 
momentum since L = —i{dg + dgr). The amplitudes ag, a±ii • • • can depend on m, but 
we make the simplifying assumption that they are independent of m. We assume that 
Gn is significant for n ~ 1, so it is sufficient to consider and a±i. Symmetry under 
9^—9 then implies that ai = a_i, and we have simply 

'ijj(9) 'ijj'{9, 9') = [oo + 2ai cos(^ - 9') + ■ ■ •] (5.57) 

We can choose oq real without loss of generality. We assume that the pion emission 
is peaked 9 — 9' , so that ai is mostly real. We now imagine that the bound 
state repeatedly emits pions, and the angular position of the pions is measured. This 
corresponds to making the transition Eq. (5.57) and then fixing 9' by picking a value 
of 9' according to the probabihty distribution 

P{9') = jd9\ilj'{9,9')\\ (5.58) 

Picking 9' in this way then gives a new wavefunction that depends only on 9, which 
can then undergo further transitions. 

This simple toy model captures the basic quantum kinematics of the problem we 
care about. For example, we can easily see how repeated transitions of the form 
Eq. (5.57) make the angular position more well-determined. The value of 9' is cor- 
related with the direction of 9, so this tends to make the peak more pronounced. 
The quantity of most interest to us is the probability to find the bound state in an 
m = state after transtions. This probability is expected to decrease as ~ 1/ V^, 
since each transition changes the maximum angular momentum by ±1. This is born 
out by Monte Carlo simulation of this model (See Fig. 8). Although this toy model 
is a drastic simplification of the system of interest, it illustrates that the expected 
behavior does arise from quantum mechanics. We therefore expect the same behavior 
in the realistic system. 
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Fig. 8. Monte Carlo simulation of toy model of angular decoherence. 
The plot shows the probability to find the system in the m = an- 
gular momentum state after interactions. The curve is the fit to 
constant / \fN . 



5.6 Hadronic Fireballs? 

The arguments above suggest that a significant fraction of colored quirk pairs lose 
most of their energy to emission of QCD hadrons. This requires that the quirks 
do not annihilate for a number of crossings of order rriQ /Aqcd ~ 10^. The non- 
perturbative QCD interactions remain effective up to very large angular momenta, 
of order £max ~ ""^q/Aqcd ~ 10^, which takes would take of order 10^ crossing 
times to reach according to the random-walk picture. In the meantime, each non- 
perturbative QCD interaction results in the emission of one (or several) hadrons with 
total energy ~ Aqcd ~ GeV. This means that the kinetic energy of the bound state 
[K ~ rriQ ~ TeV) is rapidly converted to ~ 10^ hadrons with energy ~ GeV each: a 
hadronic "fireball." 

We can obtain a simple estimate of the fraction of events of this type by assuming 
that the quirk survival probability at the n**^ crossing is (1 — P/ y/n), where P ~ 1/360 
is the s-wave annihilation probability. The probability to survive for 10^ crossings is 
then approximately 85%. 
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When the quirks finally annihilate, they are essentially at rest in their center of 
mass frame, so the annihilation products appear as a narrow resonance with mass 
2mQ. The intrinsic width will be due to the fact that the final annihilation will take 
place from a distribution of low-lying Coulombic "quirkonium" states. The width will 
therefore be of order the spacing of low- lying Coulombic energy levels, given by 

AB~afeW^~3Gev(l„^^j{^)". (5.69) 

This also sets the scale for the energy emission during the final stages of the decay, 
which wc see is only slightly larger than the QCD scale. 

The time for this process is set by the classical crossing time and the number of 
interactions required to lose the kinetic energy: 

Aqcd a2 Vmcv; Vtcv; ^ ' 

We see that the decay may have a displaced vertex for smaller values of A. 

The dominant decay will be to two jets, which may be a difficult signal due to 
large backgrounds. The decay to leptons or photons has a suppressed branching 
ratio, but offers a cleaner signal that may be easier to look for. If energy loss due to 
QCD interactions is efficient, the final hard annihilation of the quirks will be from a 
Coulomb-like state that is color and infracolor singlet. This means that there are fewer 
colored channels compared to the excited annihilation computed in Subsection 5.4. 
Assuming s-wave annihilation we find 



T{QQ jets) a 



(5.61) 



T(QQ 77) e^a^ „ 

LIA = 36^— ~ 4 X 10-^ 5.62 

r(gg^jets) al ^ > 

for eq — \. This looks very promising. The decay to is absent, although there is 
a suppressed decay mode gg^. 

Since the quirks lose all their kinetic energy before decaying in these events, the 
decay products will have an invariant mass very close to 2mQ. The intrinsic width 
will be due to the fact that the final annihilation will take place from a distribution 
of low-lying Coulombic states. These have very small energy differences of order AE 
(see Eq. (5.59)), so the intrinsic width of the resonance is very small. 

Can we hope to see the hadronic fireballs associated with these decays? Most 
of the hadrons are expected to be pions. Muons from charged pion decays will be 
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Fig. 9. Schematic depiction of hadronic fireball and hard annihilation 
into muons. Note that the the asymmetry of the muons and the fireball 
are in the same direction. 

difficult to detect because they are highly curved in the magnetic field of the detector. 
Neutral pions decay to photons, which may be more promising to detect. The angular 
distribution of the fireball may aid in distinguishing it from background. Due to the 
angular decoherence, the quirk pair acquires an angular position in the center of mass 
frame. We expect that hadron emission is peaked in the direction of the quirk motion, 
resulting in a doubly-peaked pattern in the center of mass frame. Note that the quirk 
annihilation is dominantly s-wave, and so the direction of the annihilation products 
is not correlated with the direction of the original quirk motion. This means that the 
fireball generally does not line up with the annihilation products. Furthermore, the 
longitudinal boost of the center of mass system will push both the fireball and the 
hard annihilation products in the same direction. This is illustrated in Fig. 9. 

Although we expect that energy loss due to QCD brown muck is efficient, a sig- 
nificant fraction of quirks annihilate after only a few crossings (see Eq. (5.37)). The 
branching ratio for these annihilations into leptons or photons are much smaller than 
the decays above (see Eqs. (5.40) and (5.41)), but may be worth searching for. The 
width of this enhancement is of order itlq, and the shape is determined from the 2- 
particle invariant mass distribution of the produced quirks. This gives an additional 
handle on these events. 

5. 7 Non-perturbative Infracolor Interactions 

We now consider non-perturbative infracolor interactions of the quirks. There are 
many analogies with the non-perturbative QCD interactions of colored quirks dis- 
cussed in the previous subsection, so our discussion will be brief and highhght the 
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important differences. 

The infracolor "brown muck" has a geometrical cross section for interaction, so we 
also expect ~ 1 interaction per classical crossing time. As argued in Subsection 3.3, 
radiation of infracolor glueballs takes place only while the quirk separation is less 
than or of order A~^. The non-perturbative infracolor interactions will therefore give 
rise to the emission of only ~ 1 infracolor gluons with total energy ~ A. 

One important difference with the QCD case is that the infracolor hadrons gen- 
erally do not interact after they are emitted, and therefore their angular position is 
probably not "measured" on time scales relevant for colliders. The cross section for 
an infracolor glueball with energy ~ A to scatter e.g. via — > ^g is of order 

1 4.10-.V„,(A_y(m)-', (5.63) 
167rm^ VGeVy VTeVy ' ^ ' 

where cjw ~ A^/lGvrM^ is a typical weak cross section. However, even if we assume 
that quantum coherence is maintained between the angular wavefunction and the 
wavefunction of the emitted infracolor hadrons, we still expect the probability to find 
the quirks in a low partial wave after N interactions to go like 1 / -\/]V, since the quirk 
wavefunction is "random-walking" away from low partial waves with each interaction. 
We therefore expect these interactions to suppress annihilation similarly to the QCD 
case. 

Another potentially important difference from the QCD case is the fact that the 
glueball mass is of order the strong interaction scale, so it is possible that glueball 
emission is kinematically suppressed. For example, lattice simulations of SU (3) gauge 
theory indicate that the mass of the 0"*"^ glueball is 3.6 times heavier than the square 
root of the string tension [18]. Although there is no parametric suppression, one 
should keep in mind the possibility that there is some kinematic suppression of glue- 
ball production. The amplitude to emit a hard infracolor gluon is shown in Fig. 10. 
The amphtude has one off-shell gluon and one off-shell heavy quirk line, and therefore 
is suppressed by 1/g^ where q is the hard momentum transfer, so the cross section is 
down by at large q. If this behavior sets in already at the glueball mass, we can 
imagine a suppression of order (|)^ ~ 10~^ in the cross section. 

We therefore consider the two extreme scenarios: one where there is no suppres- 
sion, and one where non-pcrturbativc infracolor interactions are effectively absent. In 
the first case, there can be significant energy loss to infracolor gluons for sufficiently 
large A, which gives rise to unobservable missing energy. In the second case, other 
mechanisms of energy loss {e.g. radiation) may be important. 
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Fig. 10. Diagram contributing to infracolor energy loss at high mo- 
mentum transfer. Hard internal lines are shaded. 



5.8 Magnetic Field 

Another effect that can be important in preventing annihilation is the magnetic field 
in the detector, of order Tesla at the Tevatron and LHC The quirk center of mass 
frame is boosted relative to the lab frame, so there will be an electric field in this frame. 
This electric field will typically have a component perpendicular to the direction of 
the quirk motion, which will give rise to a repulsive force between oppositely charged 
quirks. This in turn will generate a classical separation after one oscillation of order 

Ar~aT^~^(^)'. (5.64) 



where Vcm is the velocity of the center of mass frame. Demanding that this is larger 
than uiq^ (and assuming Vcm ~ 1) gives 

Since magnetic fields at Tevatron and LHC colliders are of order Tesla, we expect 
that this will prevent re-annihilation for A < MeV. 

This mechanism will be ineffective for special kinematical configurations. The 
induced electric field vanishes if the center of mass of the system is along the magnetic 
field, that is the beam direction. The induced electric field does not give rise to a 
transverse separation between the quirks if their motion in the center of mass frame 
is along the magnetic field. It is straightforward to check that the magnetic field in 
the center of mass frame does not cause a transverse quirk separation to leading order 
in the magnetic field. Therefore, if the magnetic field is the only effect preventing 
annihilation, there may be some events in corners of kinematic phase space that 
annihilate. 



36 



5.9 Electromagnetic Radiation 



The rate of electromagnetic radiation can be estimated from the Larmor formula 

E ~ a{df, (5.66) 

where d is the dipole moment of the charge distribution (with the charge factored 
out). A perfect s-wave has = 0, but even if the total angular momentum is f ~ 1 
the dipole moment will be of order Tmax- The energy radiated in a crossing time T is 
therefore 

AE^ETr^a ~ |. (5.67) 

Since the typical photon energy radiated is of order i?^ ~ 1/^) ^^is means that there 
are of order a photons emitted in each classical crossing. The number of crossings 
required to lose energy of order ttiq to electromagnetic radiation is therefore of order 

^ E T aA2 {GeVj \TeV J ^ ^ 

This is much larger than the number of crossings to annihilate. Brown muck in- 
teractions can prevent annihilation, but then they will be the dominant energy loss 
mechanism since E^ <^ A. We conclude that electromagnetic energy loss is unlikely 
to be important. 



6 Mesoscopic Strings 

We now consider the case where the strings are too small to be resolved in a de- 
tector (roughly L < mm), but are large compared to atomic scales (L > A). This 
corresponds to roughly 

10 keV < A < MeV (6.1) 

for niQ ~ TeV. In this case, the quirk-antiquirk pair will appear as a single particle 
in the detector. 

For mesoscopic strings, we can no longer take for granted that matter interactions 
will randomize the angular momentum and prevent the quirks from annihilating. In 
other words, we need to know whether the bound state lives long enough to appear 
in the detector. The interaction region has an inner radius of order cm with very 
high vacuum, and matter interactions are not important there. We must therefore 
consider other mechanisms to prevent annihilation. 
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An important effect in preventing anniliiiation is tfie magnetic field. In tlie quirk 
center of mass frame, there will be an electric field with a component perpendicular 
to the direction of quirk motion that gives rise to a separation of classical quirk 
trajectories after one oscillation of order 



ruQ V AV 



(6.2) 



where Vcm is the velocity of the center of mass frame. Demanding that this is larger 
than ttiq^ (and assuming Vcm ^ 1) gives 

Since magnetic fields at Tevatron and LHC coUiders are of order Tesla, we expect 
that this will prevent re-annihilation for A < MeV. 

While the quirk pair is inside the beam pipe, the only efficient mechanism for 
energy loss and change of angular momentum is the brown-muck interactions dis- 
cussed in the previous section. For colored quirks, these lead to a decay length (see 
Eq. (5.60)) 



A 

CT ~ cm ( — 
VlOl 

while for uncolored quirks 



100 keV; VTeV; 



(^y (6.4) 



We see that these decays can allow the quirk bound state to survive for distances of 
order cm. As discussed in the previous section, the efficiency of this mechanism of 
energy loss is uncertain, particularly for the infracolor energy loss. The decay lengths 
may therefore be significantly longer than these estimates. 

Once the bound state reaches the beam pipe, matter interactions are efficient at 
randomizing the angular momentum and preventing annihilation. For example, a 
single collision with an electron transfers momentum of order mg, which changes the 
angular momentum by 

For the remainder of this section we will assume that the bound state appears as a 
stable particle in the detector. In order to see the bound state, it must be produced in 
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association with a hard jet or photon so that the bound state is off the beam axis. If 
the bound state has a net electromagnetic charge, it will leave a track in the detector. 
The signal is then a single heavy stable particle recoihng against a hard jet or photon. 

The most interesting aspect of these events is the fact that the mass of the bound 
state is the invariant mass of the quirk- antiquirk pair. This has a broad distribution, 
so the mass of the bound state differs by order 1 event by event. The mass of a 
heavy stable charged particle can be measured event by event by a combination of its 
bending in a magnetic field and time of flight. This has been studied at LHC [19], with 
the conclusion that the mass can be determined at the few percent level for events 
with 0.6 < /3 < 0.8. Observation of stable particles with the mass spectrum given 
by the 2-particle invariant mass spectrum would be essentially a direct observation 
of strings. 

7 Microscopic Strings 

We now consider the signals for microscopic strings, roughly L < A, corresponding 
roughly to 

MeV < A ^ mg/few. (7.1) 

As we have seen above, in this regime interactions with matter and the magnetic field 
of the detector do not prevent the quirks from annihilating. These signals have been 
largely discussed in Section 5, so our discussion here is mainly a summary of this 
discussion. 

7.1 Colored Quirks 

We begin with colored quirks, which are the ones most copiously produced at a hadron 
collider. Most of the quirks produced above threshold will undergo hard annihilation 
without significant energy loss. Colored quirks will annihilate dominantly into jets, 
but may have branching fractions into leptons or photons at the percent level (see 
Subsection 5.4). These events will have a broad distribution essentially given by the 
perturbative 2-particle invariant mass spectrum of the quirks. 

A significant fraction (a few percent) of colored quirk pairs will lose most of their 
kinetic energy energy because of interactions of the non-perturbative QCD and/or 
infracolor interactions. The condition for non-perturbative QCD interactions to dom- 
inate is naively Aqcd > A, but there is significant uncertaintly in the efficiency of the 
non-perturbative infracolor energy loss. If the QCD interactions dominate, an energy 
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of order 2mQ will be radiated as light QCD hadrons (mainly pions) each with energy 
of order GeV: a hadronic fireball. 

The total invariant mass distribution of the quirk decay products is therefore a 
broad distribution approximating the 2-particle invariant mass distribution with a 
narrow peak superimposed. This distribution can be found in jets, but also (with 
reduced rate) in lepton or photon pairs. Detailed study of this signal would be very 
interesting. 

7.2 Uncolored Quirks 

Uncolored electromagnetically charged quirks will annihilate dominantly into infra- 
color glueballs for the range of A of interest. Infracolor glueballs are unobservable 
unless A > 50 GeV (see Eq. (2.3)), in which case they decay inside the detector. How- 
ever, there will be branching ratio typically of order 10% for annihilation into visible 

states (sec Eqs. (5.38) and (5.39)). We can have annihilation to photon pairs, or 
through a virtual photon, Z, or W . The final state will therefore include a significant 
fraction of photon and lepton pairs, which arc readily observable. 

Non-perturbative infracolor interactions will tend to bring the quirks to rest before 
they annihilate. If these are fully efficient, they will radiate of order one infracolor 
glueball with energy of order A and total angular momentum of order 1 once per 
classical crossing time. However, glueball masses are themselves of order A, so it 
is possible that there is a kinematic suppression of this process. Given our lack of 
understanding of this non-perturbative dynamics, it makes sense to consider both the 
case where these interactions are efficient and inefficient. 

Hard annihilation of excited uncolored quirks requires of order 10^ crossings (see 
Eq. (5.38)), giving the non-perturbative infracolor interactions plenty of time to ran- 
domize the angular momentum. If these interactions are efficient, wc therefore expect 
the majority of these annihilations to take place with the quirks at rest, leading to a 
very narrow resonance. If these interactions are inefficient, the resonance will be broad 
with a narrow peak superimposed from the quirks that are produced near threshold. 

8 Conclusions 

We have seen that massive particles charged under an unbroken non-abelian gauge 
group give rise to spectacular phenomenology at colliders. These signals are suffi- 
ciently exotic that they will almost certainly be missed unless they are searched for. 
Given the simple nature of these models, it is worthwhile to put some effort in this 
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direction. The next step will be to produce event generators for this exotic physics 
that can be used to develop concrete search strategies. Cosmological aspects of these 
models will also be addressed in a future publication. 
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